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Ultrathin manganite films are widely used as active electrodes in organic spintronic devices. In this study, a 
scanning tunnelling microscopy (STM) investigation with atomic resolution revealed previously unknown 
surface features consisting of small non- stoichiometric islands. Based upon this evidence, a new mechanism 
for the growth of these complex materials is proposed. It is suggested that the non-stoichiometric islands 
result from nucleation centres that are below the critical threshold size required for stoichiometric 
crystalline growth. These islands represent a kinetic intermediate of single-layer growth regardless of the 
film thickness, and should be considered and possibly controlled in manganite thin-film applications. 

Complex magnetic oxides display intriguing bulk and surface properties that make these materials very 
promising candidates for spintronics and other magnetic applications 1 . However, controlling material 
stoichiometry and homogeneity is extremely difficult, especially at the nanoscale. Therefore, the char- 
acterisation of various nanoscale features that may affect certain material parameters is important for enabling the 
manufacture of robust and reproducible devices. 

Scanning tunnelling microscopy (STM) and spectroscopy (STS) are powerful tools for characterising the 
nanoscale properties of surface electronic states, acquiring local information without significant lateral 
averaging 2,3 . 

Achieving an atomic-scale resolution in the STM imaging of manganites is an arduous task due to both the 
non-ideal atomically smooth surface topology and the high charge delocalisation 4 " 6 . The information obtained to 
date is sparse and does not provide a general or quantitative understanding of manganite surface features. To try 
to overcome these deficits, we investigated ultrathin films of the prototypical manganite La 0 .7Sr 0 .3MnO 3 (LSMO). 
Our investigation of this material was particularly motivated by its relevance to organic spintronics devices 7 " 9 . 
Thicknesses of just 2-3 nm ensure sufficient damping of the system metallicity and simultaneously decrease the 
number of statistical defects, which increases with thickness. Investigating such systems allowed us to unam- 
biguously image the atomic structures of ultrathin films of LSMO. Our results provide new insights into the 
surface properties of manganites as well as suggesting a new mechanism for the film growth of such complex 
materials. We detected a minority phase dispersed on a stoichiometric crystalline matrix; this phase consisted of 
small (less than tens of nm 2 ) non-stoichiometric islands that represent a kinetic intermediate of single-layer 
growth. We propose a phenomenological mechanism for the formation of such defects suggesting that they are 
characteristic features of all manganite films, although the detection of these defects is mainly obstructed at 
greater thicknesses. 



Results 

Fig. 1 summarises certain relevant properties of manganite ultrathin films. Resistance (R) along with magne- 
toresistance (MR) characterisations provide a comprehensive description of the electric and magnetic properties 
of these materials and reveal the characteristic metal-insulator transition (MIT). The R vs. temperature (T) 
dependence and low- field MR of the 4 samples described in Table 1 are shown in Fig. 1(a). The 15-nm-thick 
sample had no transition to a semiconducting-like behaviour in the R(T) curve within the temperature range 
investigated, and the Curie temperature of this sample as determined from the MR(T) curve was nearly 320 K. 
Such behaviour is typical of LSMO films of this thickness, making this sample a benchmark 10 . The MIT tem- 
perature (T MI ) and the Curie temperature (T c ) shifted to approximately 300 K for the 6-nm-thick film and 280 K 
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Figure 1 | (a) Resistance (R) (black curves) and magnetoresistance (MR) (red curves) versus temperature (T) for four samples of different thicknesses 
across the metal-insulator transition (MIT): 2.2 nm, 3.7 nm, 6 nm and 15 nm. (b) 1000 X 1000 nm 2 topographic image of the 2.2-nm-thick sample 
(Vbias = 1-2 V, I = 0.06 nA). (c-f) 40 X 40 nm 2 micrographs of the four sample thicknesses, respectively c) 2.2 nm, d) 3.7 nm, e) 6 nm and f) 15 nm. 
(Vbias = 1.2-0.6 V, I = 0.06-0.12 nA, depending on the sample thickness and local tunnelling conditions), (g) Room temperature differential 
tunnel conductivity curves measured on the three indicated samples. 



for the 3.7-nm sample. Finally, the 2.2-nm-thick sample showed 
neither a metallic phase nor MR. This reduction in magnetism and 
metallicity with decreasing film thickness is qualitatively and quan- 
titatively similar to results reported recently for manganite thin films 
of the highest quality, which exhibited thickness variations of 
approximately 1-2 nm between different groups and deposition 
methods 10 " 14 . 

The STM characterisations performed at room temperature (RT) 
revealed no large surface defects (such as outgrowths) for areas as 
large as 2 X 2 (am 2 . In Fig. 1(b), a 1000 X 1000 nm 2 topographic 
region of the thinnest sample displays growth terraces 100-200 nm 
wide characterised by approximately 0.4-nm height steps (approxi- 
mately one LSMO unit cell (u.c.)), in good agreement with previous 
STM reports and the layer-by-layer growth mechanism expected for 
LSMO films with small thicknesses 1215 . At this scale, the measured 
root mean square (RMS) roughness was between 0.5 and 1.0 u.c. 
and gradually increased with increasing film thickness up to 2.0- 
2.5 u.c. for the 15-nm-thick sample. The evolution of the surface 
morphology (Fig. l(c-f)) for the investigated thicknesses can be 
described in relation to the frequently observed grain-like struc- 
tures 616 , in which the topological complexity increased along with 
the thickness. 



The surface tunnelling spectroscopy is shown in Fig 1(g), where 
the dl/dV curves for three different thicknesses are plotted. All the 
spectra were averaged and acquired at the same tunnel resistance (see 
the Methods section). The curves illustrate the evolution of the sam- 
ple from a paramagnetic insulating (PI) to a ferromagnetic metallic 
(FM) phase, which is accompanied by a consistent increase in the 
metallicity, in agreement with the R(T) curves in Fig. 1(a). We also 
observed no significant correlation between the topographic and 
spectroscopic features of the films at room temperature, consistent 
with previous reports on LSMO 17 " 19 . Finally, the correspondence 
between the STS and the transport measurements confirmed that 
the physical properties of the local microscopic probes directly influ- 
enced the macroscopic device parameters. 

Further investigation of the grain-like structures revealed interest- 
ing details, especially for the two samples located in the insulating 
phase before the MIT. Fig. 2(a) displays a 20 X 20 nm 2 STM micro- 
graph of the 2.2-nm-thick sample, in which regions spatially con- 
fined within one or a few tens of nm 2 stack approximately 0.2 nm 
(one-half of the u.c.) above the surrounding "background", as shown 
by the line profile in Fig. 2(b). Such regions were characteristic of the 
entire LSMO surface, and the typical geometry of these areas is 
depicted in Fig. 2(c), wherein the 0.2-nm-step height above the grain 



Table 1 | Films thickness vs transport properties and STM findings 
Nominal Thickness [nm] Measured Thickness [nm] 



M-l Transition 



MR @ RT 



Atomic Resolution @ RT 



2 

4 

6 

13 

25 

70 



(2.2 ± 0.1) 
(3.7 ±0.1) 
(6± l) b 
(15 ±2) b 
(20 ± 2) b 
(75 ± 5) b 



NO 
<RT 
RT 
>RT 
>RT 
>RT 



NO 
NO 
YES 
YES 
YES 
YES 



YES 
YES 
NO 
NO 
NO 
NO 



a by direct XRR measurement; 

b by profilometer-based calibration. 
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Figure 2 | (a) 20 X 20 nm 2 STM micrograph of the 2.2-nm-thick sample showing a detail of the grain-like surface (Vbi as = 1.2 V, I = 0.06 nA). 
(b) Line profile over the grain surface, indicating a step height of 0.2 nm. (c) Magnification of the image in (a), where the red to blue transition 
corresponds to a step height of 0.2 nm along the z-scale, revealing a defect island. 



is indicated by different colours. We observed similar half- cell ter- 
minated regions on both the 2.2-nm and 3.7-nm samples, while the 
observation of such regions on the 6-, 15- and 75-nm-thick samples 
was prevented by slight roughness increases. 

Further magnification enabled the observation of the half-cell 
stacks at an atomic resolution, the primary result of this study. 
Fig. 3(a-b) shows atomically resolved regions of the 2.2-nm- and 
3.7-nm-thick samples. Fig. 3(a) covers an area of a few tens of nm 2 
and reveals a region of ordered atoms characterised by sizeable angu- 
lar distortion with respect to the expected square perovskite-like cell. 
A number of similar atomic patterns were detected on the same 
sample within distances of less than hundreds of nm. Among these, 
the atomically resolved regions characterised by larger areas resulted 
in a smaller distortion (Fig. 3(b)). This trend (the larger the area, the 
smaller the distortion) was confirmed on both of the thin samples. 
An analysis of this phenomenon is presented in Fig. 3(c); the inset 
indicates the method used to estimate the characteristic angle of the 
pseudo-cells. It is evident that the distortions (and hence the devi- 
ation from the standard structure) tend to disappear as the size of the 
island increases. This result is opposite to the behaviour expected 
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Island Area [nm 2 ] 

Figure 3 | (a) 5 X 10 nm 2 STM micrograph with a single defect island 
featuring an angular deviation from the expected perovskite-like square 
shape (V bias = 1.2 V, I = 0.06 nA). (b) 11 X 11 nm 2 image featuring a 
square LSMO cell (V bias = 1.2 V, I = 0.06 nA). (c) Plot of the angle vs. 
island area for the two thinnest samples; the inset indicates how we 
estimated the angle starting from a given image (the sample image is an 
offline zoom of an atomically resolved image of the 3.7-nm sample). 



from surface reconstruction effects, wherein larger groups of atoms 
result in stronger stabilisation of the distorted reconstructed sur- 
face 20 . Moreover, the distortion magnitude observed in our films 
was thickness dependent. Fig. 3(c) indicates that the changes in the 
angle of distortion as a function of island area are clearly different for 
the 2.2-nm- and 3.7-nm-thick samples. 

Discussion 

The data presented above reveal previously unknown and unpre- 
dicted features of the growth of LSMO films. In discussing these 
results, we should stress that LSMO represents a stoichiometrically 
complex compound. Considering the La/Sr molar ratio in the com- 
position of a La(0.7)Sr(0.3)MnO 3 film, it is evident that a minimum 
of 7 + 3 = 10 u.c. along two perpendicular axes in the (001) plane are 
necessary to achieve the stoichiometric composition for a layer 
grown on top of the already established stoichiometric LSMO under- 
layer. This defines an area of at least 100 u.c. or approximately 
15 nm 2 for a cell parameter of a = 0.3873 nm (see Methods). 

Furthermore, the La/Sr ratio and distribution typically deviate 
from the perfect value (7/3) due to composition fluctuations in the 
atomic flux arising from the target and due to the different La and Sr 
adatom mobilities 21 . Although such deviations have never been 
accurately quantified for thin manganite films, it is evident that the 
proposed estimation of "7 + 3 u.c." fails even for weak deviations. 
Indeed, a composition of La(0.71)Sr(0.29)MnO 3 would require a set 
of 71 +29 u.c. to guarantee the stoichiometry. Therefore, the size of 
the island necessary to establish a perovskite structure with the 
required composition/stoichiometry would be increased by one 
order of magnitude. 

These considerations led us to propose a critical island size. The 
crude estimation proposed above suggests an upper level boundary, 
such that a critical surface area S c < 150 nm 2 (1000 cells) should be 
able to stabilise the perovskite phase with the stoichiometry corres- 
ponding to the previously formed layers. This order-of-magnitude 
estimation agrees well with Fig. 3(c), wherein the distorted phase is 
observed in areas less than 120 nm 2 . 

Within this framework, adatom islands exceeding the S c evolve by 
establishing a stoichiometric perovskite structure, while smaller 
islands (S < S c ) are expected to develop temporary 2D metastable 
structures with all atoms, i.e., La, Sr, Mn and O, lying in the same 
plane and ordered in a distorted lattice, minimising the system 
energy with respect to the random orientation. We anticipate that 
this behaviour will be characteristic of any multi- composite material 
grown via Frank- van der Merwe 2D growth 22 , whereas for a step-flow 
growth, the stoichiometric defects would be confined near the frozen 
step edge. 
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This model provides a straightforward physical explanation of the 
detected atomically resolved regions. The interruption of the film 
deposition and the subsequent freezing of the structure will obviously 
result in a random number of adatoms that are unable to produce 
fully complete layers. We therefore propose that the regions observed 
at an atomic resolution correspond to defects formed by islands with 
areas that are below the critical area (S < S c ). 

The role of defects in the atomic resolution of manganites has 
already been widely emphasized: in manganites, despite the low den- 
sity of charge carriers (~10 22 cm" 3 ) compared to metals, charge 
screening is a strong limiting factor in the maximum lateral resolu- 
tion achievable by STM 4-6 . Static distortions identified as trapped 
polarons have been proposed as a necessary requisite to achieve 
atomic resolution in layered manganite La2-2jcSri+2jcMn 2 0 7 4 and 
in Lao.75Ca 0 .25Mn0 3 5 thin films. 

Moreover, atomically resolved images have been obtained only 
from the PI phase of manganites 4 ' 23 " 25 , similar to our set of samples 
(see Table 1). The semiconducting behaviour of ultrathin LSMO 
films below the MIT favours the resolution of individual atoms by 
STM because this behaviour is related to the presence of localised 
states on the surface that enhance spatial variations in the tunnelling 
current between the surface and the tip. 

The thickness dependence shown in Fig. 3(c) and the lack of 
atomically resolved islands in the thicker samples can be explained 
based upon the proposed growth model. In the 2D-island growth 
mode, the film roughness is expected to increase with increasing film 
thickness due to residual non-coalesced small islands 22 . The higher 
the density of the surface defects, the higher the density of the nuc- 
leation centres that lead to small-area islands 22 . We expect that 6- and 
15-nm-thick samples would feature a large number of a few mid- 
sized non-stoichiometric islands that cannot be resolved internally 
via STM due to both the higher system metallicity and the higher 
surface roughness. Indeed, for the thicker samples, no atomically 
resolved regions were detected, although 0.2 nm (one-half u.c.) steps 
were widely present on the film surfaces. 

Increasing the film thickness augments the metallicity of the man- 
ganite films, enhancing the screening effects in the STM visualisation 
(see above). This further decreases the ability of STM to resolve 
images with an atomic scale resolution. Our observation sheds new 
light on the mechanism of perovskite film growth by examining the 
monolayer growth process. (The cell-by-cell growth models based on 
RHEED oscillations are mere approximations). Fig. 4 illustrates the 
proposed model of the distorted, atomically resolved islands. The 
observed 2D non-stoichiometric islands represent strong potential 
defects for both charge and spin injection from the manganite elec- 
trodes. Understanding and controlling these islands is therefore of 
great importance. A deeper understanding of these effects will 
require the use of complementary surface techniques and consid- 
erable modelling efforts. However, we can rule out some possible 
competing mechanisms like surface reconstruction (see above) or 
the effects caused by the topological peculiarity of the surface, such 
as ion-reduced coordination, not able to generate half u.c. steps. 

In conclusion, ultrathin LSMO films grown epitaxially on SrTi0 3 
substrates were visualised with unprecedented size resolution. We 
emphasise the role played by inhomogeneities in the uppermost layer 
formation. Our findings revealed a separate surface phase associated 
with non-stoichiometrical defects in the uppermost film layer and 
promote significant revisions to the LSMO film growth mechanism. 



Methods 

Film deposition. High quality colossal magnetoresistance manganite LSMO films 
were grown on cubic (a = 3. 905 A) SrTi0 3 (STO) substrates using the channel spark 
ablation (CSA) technique 26 " 28 . Bulk LSMO is rhombohedral at RT with a pseudocubic 
parameter a = 3. 873 A and a distorted pseudo-cubic angle a = 90.26 029 . The quality of 
the substrates was verified with atomic force microscopy (AFM) measurements, 
which revealed terraces of ~200 nm size, with step heights of approximately 0.4 nm, 
indicating a unique surface termination. 



100 




Figure 4 | Schematisation of the "freezing" adatoms - forming a defect 
island - upon the "perfect" and mono-terminated LSMO surface. 

STM measurements. After the film deposition, the samples were transferred into the 
Omicron STM chamber under ultra high vacuum (UHV) conditions, and the surface 
was refreshed following the annealing procedures established by photo emission 
spectroscopy (PES) investigations 21 . The microscope was tip-biased and operated 
with Ptlr tips. The images shown in Figs. 1, 2 and 3 were acquired at bias voltages of 
0.5- 1 .2 V and tunnelling currents of 0.06-0. 12 nA; the tunnelling conditions slightly 
varied depending on the sample thickness. Images were acquired over an area of 512 
X 512 pixels in topographic mode and were unfiltered but levelled by a second- order 
line. 

The dl/dV curves were acquired using the current imaging tunnelling spectroscopy 
technique, in which the tip is scanned in the constant current mode to maintain a 
constant distance from the sample. We recorded the I(V) curves on a grid of 1000 X 
1000 spectra over an area of 1000 X 1000 nm 2 . Repeated measurements indicated 
that the I(V) curves, obtained as the average of tens to hundreds of curves, were 
representative of the region considered. The dl/dV characteristics in Fig. 1(g) were 
obtained by numerical differentiation and were acquired at the same reference values, 
namely 0.8 V and 60 pA. 
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